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METHODS AND SOFTWARE FOR COSMIC RAY SCINTILLATION 

STUDIES 
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Wave Propagation, USSR Academy of Sciences, 142092 
Troitsk, Moscow Region , USSR 

Methods and programs for statistical processing of 
cosmic ray intensity measured in the world-wide network 
are proposed. 

The main instrument in scintillation studies is the 
cosmic ray spectrum constructed from the intensity onser- 
vation /I/. The problem of psctrum estimation from measur- 
ed time series has a long-standing history, but the speci- 
ficities of the phenomenon of interest do not permit the 
direct use of the classical methods. One of such specifi- 
cities is that besides the appearance of regular trends, 
which can in principle be eliminated by different filtra- 
tion methods, in the periods most interesting from the 
viewpoint of the physics of processes (for instance, 

Forbush decreases, solar flares) the statistical characte- 
ristics of the process undergo essential reconstruction, 

1, e. the process besomes nonstationary. In this case the 
concept of spectrum is not defined, and the methods based 
on Fourier transform, Blackman-Tusky method) give wrong 
results. The technique usual in this situation, i.e. discri- 
mination of quasistation ary regions, encounters the follow- 
ing diff iculties . Such regions, if any, can be very short. 
When'rbased on a small amount of datatj the methods using 
Fourier transform are known to give inaccurate results, 

in particular, they do not allow us to separate close fre- 
quencies, whereas the second specific feature is just the 
necessity to separate close frequencies, some of which can 
be associated with a certain process in solar space . For a 
better separation of near frequencies in short regions, the so- 
called autoregressive methods have lately been used for spec- 
trum estimation /2/. These methods suggest introduction of an 
additional possibility to describe the process by an auto- 
regressive model ,o 

SC. . = 5X X I + 5 Lzoj 
^ i =o t 1 t (I) 

where is the sequence of independent random quantities 

of a certain unkhown order p. This assumption is based on 
the fact that a wide class of stationary processes can be 
described by such a model. The coefficients of regression 
are in this or that way estimated tinder this assumption for 
a certai chosen order p, and then the single-va'lued ana- 
lytical estimate of the spectrum is found from these coeffi- 
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dents. Such an approach, using different algorithms for 
estimating the coefficients (of the type of Berg, Levinson- 
Durbin, Pisarenko, and other methods), is also realized in 
our methods. In some cases it gives satisfactory results. 
But it cannot be applied to essentially nonstationary pro- 
cesses either. ’ 

To eliminate these difficulties we have proposed the 
following two approaches. One of them suggests the descrip- 
tion of a process by an autoregressive model in which the 
coefficients vary in time 

+ (2) 

Bach coefficient is represented as a series in a given 
complete set fo functions 

= ' (3) 

with the unknown coeff icients^C^. As a set of functions 

one can choose, in particuler, a” power series 4 V r • • y. 
Then, using the method of least squares, one can calculate 
the coefficients for the chosen order p of the 

model and for the number N in the expansion (3) , The order 
p of the model and the number of terms N in the expansion 
ca be chosen, so to say, in the best possible way. 


Such an approach makes it possible to introduce the 
concept of an instantaneous spectrum for a nonstationary 
process. At each time moment (estimated from the time 
series of lenght T) the parameters correspond to_an/v . n 

autoregressive model with known constant coeff icient oJ 


oe 

l ~0 T 


£-■1 


Such a process will be thought of as stopped at a time mo- 
ment t (starting at some initial point X c , it passes 
through a stationary sequence) . This process corresponds 
to a psectrum which can be calculated analytically from 
the moment t will be called a t -instantaneous spectrum 
of the non station ary process (2) . 

Arranging the sequence of instantaneous spectra with 
respect to t , one obtains a dynamical picture of re- 
construction of a nonstationary process. 

The other approach to investigation of a nonstationa- 
ry process is based on the on-iine method for estimating 
the parameters of the model (2), i.e. on the stochastic 
approximation method. Let the coeff icient s{ oC[ (i*)Jbe time- 
inde pendent (model (I)) , then the algorithm of stochastic 
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approximation for estimating ia v \ for a chosen order p 
of the model has the form (a is the vector of parameter 
estimation, vp - ( x, , • • • , x ) is the measured time se- 
ries): * p 

flj.i = a* - Ht % (h.L - ?r c ) (4) 

Under certain conditions, if the sequence tends to zero 
with an appropriate velocity, one can show (see, for ins- 
tance, /?/) that as t co, the estimates of a* converge 
to the vector of the coefficients a of model (I). It is a 
disadvantage of the method that it needs a large sample, 
which is hardly to be provided under the described conditi- 
ons. At the same time, this difficulty can be eliminated if 
we reject the requirement concerning constancy of the coef- 
ficients and pass over to the model (2) without expanding 
a in the series (3). 

In this case, to estimate the coefficients of the model 
(2), we use the algorithm /4/, in which the sequence V , 
is chosen differently. One can show that under certain con- 
ditions on the growth of the coefficients & ' L (t) and for a 
corresponding rate of convergency ^ to a certain constant, 
the algorithm (4) will trace the trend of the parameter <Jt (t). 
Under these assumptions one can provide a large sample of 
data since the series need not necessarily be stationary, 
and the algorithm works well. The way of constructing an in- 
stantaneous spectrum is the same as described above. 

The software considered includes all the approaches 
described above, namely, direct Fourier transform and its 
modifications, autoregressive and instantaneous spectrum me- 
thods. Being used in various combinations, they prove help- 
ful in handling the time series. The methods have been veri- 
fied in special tests. Some results of investigation of scin- 
tillations using these methods are presented in /4/. 
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